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Abstract 



> 

We investigate the effect of the dimuon CP asymmetry from the B decay modes, recently 
observed at 3.2 a deviation from the Standard Model (SM) by the DO collaboration, in the 
context of SU(5) and SO(10) GUT models. We exhibit that a large amount of flavor violation 



between the second and the third generation is generated due to the large neutrino atmospheric 
mixing angle and this flavor violation can be responsible for the observed large CP asymmetry 
due to the presence of new phases (not present in the CKM matrix) in the Yukawa couplings. 
We also study the implication of the parameter space in these GUT models with large CP 
violating lepton asymmetry for different phenomenologies, e.g., Br(r — \ /17), Bi(B s — > fifi) at 
the Fermilab, direct detection of dark matter (DM) in the ongoing detectors and measurement 
of muon flux from solar neutrinos at the IceCube experiment. 



1 Introduction 

Recently, the DO Collaboration has announced the appearance of a like-sign dimuon charge 
asymmetry in the semileptonic 6-hadron decays measurement: A b sl = — 0.00957±0. 00251 (stat)± 
0.00146(syst) pp. In the Standard Model (SM), the prediction for the asymmetry is A^(SM) = 
(— 2.3lo;g) x 10 -4 , and the DO experimental result differs from this by 3.2 standard deviation. 
In the absence of CP violation this quantity clearly vanishes, hence the DO result leads us to 
a new physics (NP) which induces some CP violating flavor changing interactions beyond the 
SM. 

The BB mesons created in pp collisions include both Bd(db) and B s (sb). The quantities of 
the Bd system are well measured by 5-factories, and the unitarity triangle seems to be closed: 
VudV^ b + V c dV* b + V t dV t l = 0. In that case, the asymmetry from the Bd-Bd oscillation is expected 
to be tiny and a new CP violating phase should show up in the B s system instead, namely in the 
b-s transition. The DO and CDF Collaborations have reported the existence of a CP violating 
phase, <p s , in the B s system from the B s —¥ J/ip4> decay [21 [3]. In fact, their results differ from 
the SM prediction in a direction which is consistent with the signature of the like-sign dimuon 
asymmetry reported by the DO. This provides us an encouraging guide, pointing towards a 
source for new CP violation in the b-s transitions. 

Supersymmetry (SUSY) is a very attractive candidate to build NP models. As it is well 
known, SUSY models have a natural dark matter candidate which is a neutralino as the light- 
est SUSY particle (LSP). Besides, the gauge hierarchy problem can be solved and a natural 
aspect of the theory can be developed from the weak scale to the ultra high energy scale. In 
fact, the gauge coupling constants of the Standard Model gauge symmetries can unify at a 
high scale using the renormalization group equations (RGEs) of the minimal SUSY standard 
model (MSSM). This indicates the existence of a grand unified theory (GUT) as the underlying 
principle for physics at the very high energy. The well motivated SUSY GUT models have al- 
ways been subjected to intense experimental and theoretical investigations. Identifying a GUT 
model, as currently is, will be a major focus of the upcoming experiments. 

In SUSY models, the SUSY breaking mass terms for the squarks and sleptons must be 
introduced, and they provide sources for flavor changing neutral currents (FCNCs) and CP 
violation beyond the Kobayashi-Maskawa theory. In general, soft breaking terms generate too 
large FCNCs, hence a flavor universality is often assumed in the squark and slepton mass 
matrices to avoid large FCNCs in the meson mixings and lepton flavor violations (LFV) [I]. 
The flavor universality is expected to be realized by the Planck scale physics. However, even 
if the universality is realized at a high energy scale such as the GUT scale or the Planck scale, 
non-universality in the SUSY breaking sfermion masses is still generated through the evolutions 
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of the RGEs, and this can lead to some small flavor violating transitions which could possibly 
be observed in the ongoing experiments. In some MSSM models with right-handed neutrinos, 
the induced FCNCs from the RGE effects are not large in the quark sector, while sizable effects 
can be generated in the lepton sector due to the large neutrino mixing angles [5]. Within 
GUTs, however, loop effects due to the large neutrino mixings can induce sizable FCNCs also 
in the quark sector since the GUT scale particles which connect the quark and lepton sectors can 
propagate in the loops [5J. As a result, the patterns of the induced FCNCs highly depend on the 
unification scenario and the heavy particle contents. Therefore, it is important to investigate 
FCNC effects to obtain a footprint of the GUT physics. 

If the quark-lepton unification is manifested in a GUT model, the flavor violation in the 
b-s transition can be responsible for the large atmospheric neutrino mixing [7], and thus, the 
amount of flavor violation in the b-s transition (the second and the third generation mixing), 
which is related to the B s -B s mixing and its phase, has to be related to the r — >• nq decay 
[HI El QUI El 021 02] for a given particle spectrum. The branching ratio of the r — >• /ry decay 
is being measured at the S-factory, and thus, the future results on LFV and from the ongoing 
measurement of the B s -B s mixing phase will provide important information to probe the GUT 
scale physics. 

In Refs. [HI QUI 02] , the authors studied the correlation between the branching ratio of r — > fi'-f 
and the phase of the b-s transition in the frameworks of SU(5) and SO(10) GUT models. While 
performing the analysis, it is important to pay attention to the dependence on tan/3, which is 
the ratio of the vacuum expectation values of up- and down-type Higgs bosons. In the case of 
tan/3 < 20, the gluino box diagram can dominate the SUSY contribution to the B s -B s mixing 
amplitude, while for large tan/3 ^ 30 it can be dominated by the double penguin diagram 
contribution [TJ1 [151 OS] ■ When the Dirac neutrino Yukawa coupling is the origin of the FCNC 
(we call this case as minimal type of SU(5)), the r — > \x^ constraint gives a strong bound on the 
phase of the B s -B s mixing for smaller tan/3. On the other hand, when the Majorana neutrino 
Yukawa coupling is the source for the FCNC, both left- and right-handed squark mass matrices 
can have off-diagonal elements (we call this case as minimal type of SO(10)), the gluino box 
contribution is enhanced and larger B s -B s phase is possible compared to the SU(5) case. The 
double penguin contribution is proportional to tan 4 /3, while the Br(r — > fi'-y) is proportional to 
tan 2 /3. As a result, for both SU(5) and SO(10) cases, a large phase of b-s transition is allowed. 
In that case, however, Br(B s —¥ h/j,) constraint is more important since it is proportional to 
tan 6 /3 [17]. In other words, when the phase of the b-s transition is large due to the double 
penguin contribution, Bt(B s — > /i/x) has to be also enhanced. In fact, in [12] it was shown that 
a lower bound from Br(B s —¥ fi/i) is obtained in the case of SU(5) GUT. 



In Ref . [18] , we also investigated the implication on the dark matter detection from the large 
B s -B s mixing. Assuming that the neutralino LSP is the dark matter candidate, the SUSY 
parameters are restricted by the relic density constraint. It was shown that the funnel region, 
in which the relic density constraint is satisfied through annihilation near heavy Higgs pole, is 
favored by the flavor solution. Moreover, there are some correlation between flavor changing 
processes and the neutralino direct detection cross section through the dependency on the 
CP-odd Higgs mass, m^. 

In this paper, we will sort out the GUT models, where the FCNC is due to the atmospheric 
neutrino mixing, to obtain a large CP asymmetry of the B decays. This investigation is 
important if the reported size of the like-sign dimuon charge asymmetry persists in the future 
with a smaller error. We show that it prefers the SO(10) type boundary condition where 
Majorana neutrino couplings induce the FCNC and both left- and right-handed squark mass 
matrices have off-diagonal elements. Especially, for the SU(5) boundary condition where the 
Dirac neutrino Yukawa coupling induces the FCNC to produce a large CP asymmetry, a large 
value of tan (3 is required and the SUSY mass spectrum is restricted. We will also study the 
implication of the dark matter direct and indirect detections from the constraints on the SUSY 
mass spectrum. 

This paper is organized as follows: in section 2 we discuss the resent results of CP violation 
in B s decays; in section 3, we discuss the sources of flavor changing neutral currents (FCNC) 
in the context of SUSY GUTs; in section 4, we show constraints arising from the experimental 
constraints on different FCNC processes; in section 5, we discuss the constraints from the dark 
matter content of the universe and predictions related to the direct and indirect detection 
experiments; and we conclude in section 6. 

2 CP violation in B s decays 

The dimuon like-sign asymmetry A h sl by DO deviates from the SM prediction by 3.2 a . The 
BB samples created at the Tevatron includes both Bd and B s , and the asymmetry can be 
written as A b sl = (0.506 ± 0.043)4 + (0.494 ± 0.043)^ pQ. The pieces of a q al {q = d, s) can be 
defined as a^ = (r q — f q )/(r q + f q ) where r q and f q are the ratios of B-B mixings: r q = P(B — > 
B)/P(B — > B) and f q = P{B — > B)/P{B — > B). Since the B d system is consistent with 
experiments, we assume that the CP asymmetry in the Bd-Bd mixing is negligible. When we 
take the experimental data on dimuon asymmetry by CDF (1.6 fb _1 ) and on "wrong-charge" 
asymmetry in the semileptonic B s decay by DO into account, the asymmetry in B s -B s mixing 
is extracted as followsll9j 

a 8 Bl = (-12.7 ±5.0) x 10~ 3 , (1) 



which deviates from the SM prediction by about 2.5 a. 

When r^ 2 <C M{* 2 (Mf 2 is the mixing amplitude and I^ 2 is the absorptive part of the mixing) 
a|j is given as [21] 
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where <p s is the argument of Y\ 2 /M[ 2 . 

In many NP models, the AB = 2 (B is beauty) Hamiltonian can be easily modified (e.g. 
see for recent works motivated by the DO results [191 20J). We parameterize the M{ 2 as 
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where C s is a real positive number. From the measurement of the mass difference, AM S = 
2|Mf 2 |, the experimental result is consistent with C s = 1. When the Ai? = 1 Hamiltonian is 
same (allowing modification at 10% level) as the SM, even in the presence of new physics, the 
phase of I^ 2 is almost same as that of SM, which is tiny ~ 0.04 (in usual phase convention 
where V c bV* s is almost real). In that case, (p s is the same as the phase (— 2/3, = — 2(/3f M + </>b s )) 
measured by B s —¥ J /ip(f) decay observation. Using the B s — > J /ip4> decay, the decay width 
difference Ar s = 2 1 Tf 2 1 cos0 s is also measured. The parameters of B s -B oscillations, /3 S and 
AT S , have been measured at the Tevatron [2], and the CDF Collaboration showed their recent 
analysis till 5.2 fb _1 of data [3]. It appears that the data statistics is very different over the 
periods (0 - 2.8 fb" 1 and 2.8 - 5.2 fb" 1 ). We will adopt their analysis for - 5.2 fb _1 . The 
CDF result on the phase 2f3 s differs from the SM prediction 2/3^ M ~ 0.04 at 1 a level (DO shows 
about 2 a deviation for the same measurement for 2.8 ftr 1 data |2J), and the signature of the 
phase is consistent with the sign required to explain the anomalous like-sign dimuon charge 
asymmetry by DO. 

The SM prediction on r* 2 /Mf 2 is given by Lenz-Nierste [22] 
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It was pointed out that theoretical prediction of the absolute value of a s sl is bounded from above 
and the bound is a little bit too small to explain the dimuon asymmetry by DO if T\ 2 is not 
modified [T91I23]. This is because the experimental measurement of AM S = 2|Mf 2 | is consistent 
with the SM prediction (which means C s ~ 1). Using the simple relation: 
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we illustrate the current situation in the Figure 1. The solid circles correspond to the 1 a 
boundaries in the case of C s = 1 by using the SM prediction by Lenz-Nierste. The colored 
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Figure 1: The experimental and theoretical regions in the a^ r AT s /AM s plane. The yellow 
region is the 1 a region of the combined data from the semi-leptonic B decays. The red and 
blue lines are 95% and 68% boundaries from the B s — > J/ip(f> decay, assuming that the phase 
of r* 2 is the same as the phase of T^ 2 sm- The solid circles are the theoretical 1 a boundaries 
using the numerical calculation in Ref.[22]. The dotted circle corresponds to the conservative 
theoretical region when one implements the remedies 2 and 3 in the text. 



(blue and red) solid curves correspond to the measurement from the B s — > J/i/xj) decay by 
CDF (0 — 5.2 fb _1 ). We assume that <fr s = —2(pB a , which means that argT^ 2 = argT^ SM . The 
horizontal colored (yellow) band is the 1 a region of a\ v As one can see, the 1 a regions do not 
match well, independent of the choice of <f> s , where tan0 s = a s sX /{AT s / AM S ). 

For the current situation, the combined analysis has a large error for a s sh the discrepancy 
is not very serious if the phase <fi s is 0(1) rad. However, the dimuon asymmetry measured by 
the DO alone has a large center value (which corresponds to a s sl = (—19.4 ± 6.1) x 10~ 3 ), and 
if the measurement of the dimuon asymmetry (or the "wrong-charge" muon asymmetry in the 
semileptonic B s decays) becomes accurate in the future keeping the large center value of af b 
one has to resolve this issue. To implement such possible future constraints, one can consider 
the following three typical remedies. 

1. Add aAB = l effective Hamiltonian to modify T s 12 or Tf 2 p4ll25| l26| 127] . 

This is a direct resolution of this issue. If in the future measurements the phase of 
B s — > J/ip(f) is really diminished, this type of resolution will be needed. (In fact, the 
recent CDF data for the period 2.8 — 5.2 fb _1 may indicate that the B s — > J/ip<f) phase is 
almost zero.) 



2. Non-perturbative effects [23| 126]. 

In this case, the numerical number in Eq.(JH) is obtained by a two parameter expansion, 
Aqcd/^6 and a s (rrib), using operator product expansion and heavy quark expansion. In 
fact, it is known that non-perturbative effects may dominate in the D°-D° meson mixings, 
and it may be true that there is a large long distance contribution in the T s 12 (e.g. the 
intermediate states include DfDj, which may lead to large non-perturbative effects). In 
the case of B s -B s mixings, each term of the next to leading order is about 30% of the 
leading order, and the expansion may be more reliable than in the case of D-D mixing. 
However, a careful treatment is needed since the series may not be converging. Actually, 
the term which has the largest uncertainty in the next to leading order calculation gives a 
negative contribution to T s l2 , and the true numerical value may be larger than in Eq.fllJ). 
In that sense, the discrepancy is not so serious unless the B s — > J/ipcj) phase will become 
tiny with a small error in future. 

3. Use the uncertainty of the Bag parameter Bb s and the decay constants /b s . 

The mixing amplitudes are proportional to BbJ'b > which has about 40% error. This 
factor is canceled in the ratio of Y\ 2 jM{ 2l and the SM prediction in Eq.(j4]) does not have 
the ambiguity. The parameter C s in Eq.Q may have the 40% error. However, since the 
ratio of the hadronic quantities for Bd and B s , related to the SU(3) flavor violation, is 
more accurate 1281 



£= ^y_gi = 1.23 ±0.04, (6) 
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the mass difference of Bd restrict the uncertainty of C s less than 10% because of the 
relation 
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Therefore, if one uses the full uncertainty of BbJb i one also has to modify \Mf 2 \ in the 
same rate as of |Mf 2 |. In general, it is possible to do that, but one should be careful about 
the argument of Mf 2 since the sin 2/3 measurement is consistent with the SM. 

In SUSY models, the argument of Mf 2 is related to the 13 mixing/23 mixing of the 
squark mass matrices, and in models where the FCNC is induced by the Dirac/Majorana 
neutrino Yukawa couplings, it is related to the size of the 13 neutrino mixing. 

In this paper, we consider the case where T\ 2 ~ T\ 2 SM and the phase (f) s comes from the 
Mf 2 phase <p s = —2(pB a , in SUSY models with i?-parity conservation. 



The dotted circle in the Figure 1 corresponds to the region when we used the 2 a error of 
the Eq.(j4]) and 40% error from the B Bs f Bs and T^ 2 = T^ SM . The absolute value of the phase 
4> s should be large ~ 50° — 70° to explain the large CP asymmetry oj. 

By definition in Eq.([I)]), we obtain 

• 2 , (af») ~ ( 1_ ^) 2 fa . 

sm <Pb s = 4C > ( 8 ) 

where A^ p = \M^ NP \ and A S S M = \M^ SM \. When C s ~ 1, we obtain 2sin0 Bs ~ A^/Af 4 . 
Therefore, the NP contribution of Mf 2 should be comparable to the SM contribution to obtain 
the large phase 4> s ~ —2<p Bg . 

In SUSY models (for an early study of the dimuon asymmetry in SUSY models, see |29j). 
we need to realize A^ p ~ A^ M in order to explain the current combined data of CP asymmetry 
in B decay. As it is mentioned already that in GUT models, the Dirac/Majorana neutrino 
Yukawa coupling can be a source for FCNC even in the quark sector. When the quark-lepton 
unification is manifested, the amount of A^ p is related to the lepton flavor violation r — > /Z7, 
and will be bounded by the constraint on Br(r — » /ry). The main purpose of this paper is to 
study how to obtain a large A^ p when there is quark-lepton unification and after satisfying all 
the other experimental constraints. 

3 FCNC sources in SUSY GUTs 

In SUSY GUT theories, it is often assumed that the SUSY breaking sfermion masses are flavor- 
universal, but the off-diagonal elements of the mass matrices are generated by the loop effects. 
The FCNC sources are the Dirac/Majorana neutrino Yukawa couplings, which are responsible 
for the large neutrino mixings [5j |6] . Since the left-handed leptons (L) and the right-handed 
down-type quarks (D c ) are unified in 5, the Dirac neutrino Yukawa couplings can be written as 
Y U ij5iNjH$, where iV c is the right-handed neutrino. The flavor non-universality of the SUSY 
breaking D c masses is generated by the colored Higgs and the iV c loop diagram [7], and the 
non-universal part of the mass matrix is <5MjL — — ^ (3rn,Q + A 2 ^)Y V Y}, ln(M*/ikfff ), where M* 
is a cut-off scale (e.g. the Planck scale), Mh c is a colored Higgs mass, m is the universal scalar 
mass and Aq is the universal scalar trilinear coupling. The left-handed Majorana neutrino 
coupling LLAl (Ax, is an SU(2)x triplet) can also provide contributions to the light neutrino 
mass (type II seesaw [30]), and can generate the FCNC in the sfermion masses when the 
fermions are unified. 

As a convention in this paper, we will call the model with the FCNC source arising from the 
Dirac neutrino Yukawa coupling as the minimal type of SU(5). In this case, the off-diagonal 



elements of 10 (Q, U c , E c ) representations are small because they originate from the CKM 
mixings. In a competitive model which we call the minimal type of SO(10), the Majorana 
couplings, which contribute to the neutrino mass, generate the off-diagonal elements for all 



sfermion species since the Majorana couplings fijLiLjAL can be unified to the fijlGi 16 j 126 
coupling [3Tj . 

We note that when the source is not specified, such as the case for Dirac Yukawa coupling, 
the off-diagonal elements (of 10 multiplets in SU(5)) can be large in general. In our convention 
of SU(5) and SO(10) models, the source of the off-diagonal elements are specified, and only 
the right-handed down-type squark mass matrix has sizable off-diagonal elements in SU(5), 
while both left- and right-handed squark mass matrices can have sizable off-diagonal elements 
inSO(10). 

The non-universal part generated from the Dirac/Majorana couplings, Y v and /, is propor- 
tional to Y V YJ and //'. In general, therefore, the squark and slepton mass matrices due to the 
loop correction can be parameterized as 

M% = mjj[l - K F U F disig(k h k 2 , 1)U F ], (9) 

where F = Q, U c , D c , L, E c . The unitary matrices Up is equal to the neutrino mixing matrix in 
a limit [32]. The quantity Kp denotes the amount of the off-diagonal elements and it depends 
on the sfermion species. In the minimal type of SU(5), since the off-diagonal elements of the 
SUSY breaking mass matrix for the left-handed lepton doublet get the correction, 8M? ~ 
— l/(87r 2 )(3mQ + Aq) J2k(Yv)ik(Y u )jk ln(M*/Mfc) where M* is a scale that the flavor universality 
is realized, k l can be written as k l ~ (Y I/ diag )^ 3 (3 + Al/ml)/(8ir 2 ) \n(M„/M R ), where M R is a 
Majorana mass of the right-handed neutrino. If the GUT threshold effects are neglected, we 
have Kg = kl = kdc 1 and U% = Ul = Ub c - In general, the fermion mass matrices arise from 
the sum of the Yukawa terms and the equality of Ub c and Ul can be completely broken when 
there are cancellations among the minimal Yukawa term and additional Yukawa terms. Here, 
we consider a model where the (near) equality between Ud° and Ul (especially for 23 mixing 
angle of them) is maintained as a "minimal type" assumption. The assumption is natural if 
there is a dominant Yukawa contribution and corrections to fit realistic masses and mixings 
are small. The unitary matrices for Q, U c , E c species are related to the CKM matrix, and can 
generate only negligible effects to the following discussion. Therefore, for the SU(5) boundary 
condition, we assume, is as follows: 

SU(5) : k l = K D o, U L = U D c, kq = «yc = Kb* = 0. (10) 

This boundary condition will be used for discussions in the following sections. In the minimal 
type of SO(10) model, all Up can have large mixings responsible for the neutrino mixings. If 
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the threshold effects are neglected, one finds k 16 ~ 15/4(/3 d 3 iag ) 2 (3 + Al/ml)/(8n 2 ) In M*/M v , 
where Mjj is a SO(10) unification scale. In general, however, the equality of all k e in the 
SO(10) boundary condition is broken by threshold effects. A detail physical interpretation of 
this parameterization is given in [12, 32J. The SO(10) boundary condition, we assume, is as 
follows: 

SO(10) : Kq = K V c = K D c = K L = K E c, Uq = U V c = U D c = U L = U E o. (11) 

When the Dirac neutrino Yukawa coupling Y v or the Major ana coupling / is hierarchical, we 
obtain k%, hi -C 1 and then the 23 element of the sfermion mass matrix is — l/2mQ/tsin2#23e ia . 
The magnitude of the FCNC between 2nd and 3rd generations is controlled by k sin 2# 23 , where 
#23 is the mixing angle in the unitary matrix. The phase parameter a also originates from the 
unitary matrix, and it will be the origin of a phase of the FCNC contribution. 

As it is mentioned that it is preferable to modify the absolute value of B^-B^ mixing am- 
plitude \Mf 2 \ without modifying its argument in order to enhance the asymmetry a s sl . The 13 
element of the sfermion mass matrix is k(— 1/2/^ sin 29\2 sin #23 + e* 5 sin #13 cos #23) [32], where 
6ij are the mixing angles and 5 is a phase in the unitary matrix. Choosing small values for the 
parameters k-i and #13, one can realize the preferred situation. 

4 Constraints from the FCNC processes in SUSY GUTs 

In the MSSM with flavor universality, the chargino box diagram dominates the SUSY contri- 
bution to M( 2 . In the general parameter space of the soft SUSY breaking terms, the gluino 
box diagram can dominate the SUSY contribution for a lower tan/3 (i.e. tan/3 < 20 — 30). The 
gluino box contribution is enhanced if both left- and right-handed down-type squark mass ma- 
trices have off-diagonal elements [33], and therefore, it is expected that the SUSY contribution 
to the B s -B s mixing amplitude is large for the SO(10) model with type II seesaw, compared to 
the minimal type of SU(5) model [10] . 

When the lepton flavor violation is correlated to the flavor violation in the right-handed 
down-type squark, the r — >• /X7 decay gives us the most important constraint to obtain the 
large B s -B s phase pTJ] [TT]. Furthermore, the squark masses are raised much more compared to 
the slepton masses due to the gaugino loop contribution since the gluino is heavier compared 
to the Bino and the Wino at low energy (assuming the gaugino mass universality at a high 
energy scale such as the GUT scale), and thus the lepton flavor violation will be more sizable 
compared to the quark flavor violation. The current experimental bound on the branching ratio 
of r — > )U7 is [34] 

Br(r -»■ /i 7 ) < 4.4 x 10~ 8 . (12) 
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Figure 2: The possible SUSY contributions are plotted when the r — > /17 bound is saturated. 
The SO(10) boundary condition can give a larger SUSY contribution than the SU(5) boundary 
condition. We choose mi/2 = 300 GeV and tan/3 = 10 for this plot. The detail to draw the 
plot is written in the text. The relation between the CP phase <\>. 
given in Eq.(jHJ). 



NP / /i SM 



-2cj) B3 and A™ /A 



is 



In order to allow for a large phase in the B s -B s mixing, a large flavor-universal scalar mass 
at the cutoff scale is preferable. The reasons are as follows. The gaugino loop effects are flavor 
invisible and they enhance the diagonal elements of the scalar mass matrices while keeping 
the off-diagonal elements unchanged. If the flavor universal scalar masses at the cutoff scale 
become larger, both Br(r — > yrf) and <J)b s are suppressed. However, Br(r — > 7x7) is much more 
suppressed compared to 4>b s for a given k sin 26*23 because the low energy slepton masses are 
sensitive to uiq while the squark masses are not so sensitive due to the gluino loop contribution 
to their masses. The large Higgsino mass, /1, is also helpful to suppress the dominant chargino 
contribution of r — > fi'y. The subdominant neutralino contribution, however, will become large 
when /i is large due to the large left-right stau mixing. 

In figure 2, we plot the magnitude of A^ p /Ag as a function of m^ (the SUSY breaking 
mass of 5 = (D C ,L) at the unification scale), when the r — > \v) bound, Eq. (fl2|) . is saturated, 
for various mass parameters in the case of tan/3 = 10. We choose the unified gaugino mass as 
mi/2 = 300 GeV, and the universal scalar trilinear coupling as A = 0. In the case of SU(5), 
the SUSY breaking mass of 10 = (Q, U c , E c ), m 10 , can be different from m^. As one can see, 
in order to achieve A^ p /A^ ~ 1, the mass parameters should be around 2 TeV. In the case 
of SO(10), the sfermion masses are unified, mo = m^ = mio, and thus we only change \x in the 
two plots in the figure. As one can see, the NP contribution in SO(10) can be much bigger than 
SU(5) case. This is the consequence of the fact that both left- and right-handed squark mass 
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matrix can have sizable off-diagonal elements from the Majorana neutrino coupling in SO(10) 
case. 

In the lower tan/3 case, however, the amount of non-universality k has to be large > 0.3 
to achieve A^ p /A^ M ~ 1, especially in SU(5). Such a sizable n is possible if A /mo is large, 
but the large k is not preferable as long as it is the RGE induced origin from the Planck scale 
and the GUT scale. Besides, the muon g — 2 anomaly |35j is also suppressed when r — > [irj 
is suppressed. This is not good since the deviation of g — 2 from the SM prediction is now 
estimated about 3.2 [3B] - A a level [37]. When tan/3 is larger (> 30 — 40), the so called double 
Higgs penguin diagram dominates the contribution rather than the box diagram, and k can be 
smaller < 0.1 to achieve A^ p /A^ ~ 1. In this case, we do not need to suppress r — >■ /ry, and 
the muon g — 2 anomaly can be explained. 

The double Higgs penguin (flavor changing neutral Higgs interaction) is generated as follows 
[17J. In SUSY models, only the holomorphic coupling is allowed for the Yukawa coupling. When 
SUSY is broken, the non-holomorphic coupling is generated by the finite corrections. For the 
down-type quark, the Yukawa coupling is 

£ eff = Y d QD c H d + Y^QD C H*. (13) 

The second term is the non-holomorphic term, and YJ23 and YJ32 are roughly proportional to 
tan/3. Since we work on the basis where the down-type quark mass matrix {M d = Y d v d + Y d v u ) 
is diagonal, the following flavor changing Higgs coupling can be obtained: 

Y>QD C H: - Y^QD c H d . (14) 

^d 

The dominant flavor changing neutral Higgs coupling is roughly obtained from the second term, 
and it is proportional to tan 2 /3. The B s -B s mixing can be generated from a double penguin 
diagram, the mixing amplitude is proportional to tan 4 /3. Since the Br(r — > ^7) is proportional 
to tan 2 /3, the double penguin contribution for large tan/3 is preferable to obtain A^ p /A^ M ~ 1 
satisfying the r — > /i7 constraint in GUT models. 

The effective flavor changing Higgs couplings are written as 

X^&Pnd^S + X^&PlcI^S , (15) 

where S° represents for the neutral Higgs fields, S = [H, h, A] , where H and h stand for heavier 
and lighter CP even neutral Higgs fields, and A is a CP odd neutral Higgs field (pseudo Higgs 
field). The couplings are 

x rl = Y 'dij^K r[sin(a-/3),cos(a-/3),-i], (16) 

V2cos/3 

x lr = Y dji^K -[sin(a-/3),cos(a-/3),z], (17) 

V2cos/3 
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Q. 




tany? 

Figure 3: The possible SUSY contributions are plotted as a function of tan /3 when the r — y /ry 
bound is saturated. The SO(10) boundary condition can give a larger SUSY contribution than 
the SU(5) boundary condition. We choose mi/2 = 300 GeV and m = 800 GeV. The detail to 
draw the plot is written in the text. 



where a is a mixing angle for h and H. The double penguin diagram including both left- and 
right-handed Higgs penguin which is proportional to the factor 



sin 2 fa 



/3) cos 2 (a 



m 



m 



2 ~T ^2 ' 



m' 



;is) 

b H " b h " b A 

and the double penguin contribution is naively proportional to Xj^X 2 ^/ 'm 2 A (HI HSl US] . On 
the other hand, the double left-handed (or double right-handed) penguin contribution oc (X^ R ) 2 
(or (Xj^) 2 ) is tiny because cos(a — /3) ~ and m A ~ m# for tan/3 3> 1 and m A > m.#. 

In the flavor universal SUSY breaking, the right-handed penguin coupling X R \ is tiny, and 
the double penguin contribution cannot be sizable even for a large tan/3. However, when the 
right-handed mixing is generated in the SUSY GUT models, the double penguin diagram can 
be sizable for large tan/3. We note that if there is a FCNC source in the right-handed squark 
mass matrix, we do not need the off-diagonal elements in the left-handed squark mass matrix 
in order to generate the sizable double penguin contribution. Therefore, even in the minimal 
type of SU(5) model, the double penguin contribution can be sizable when tan/3 is large. 

In figure 3, we plot the tan/3 dependence of A^ P /A^ M when Br(r — y fj,j) saturates the 
experimental bound for two Higgsino masses /i. We choose the unified gaugino mass as m\/2 = 
300 GeV, and the SUSY breaking scalar masses as m = m 5 = mw = 800 GeV, and A = 0. We 
assume m 2 H = m 2 H for the SUSY breaking Higgs mass at the unification scale just for reducing 
the number of parameter. As one can see, A^ p /A^ M becomes smaller for tan/3 ~ 20. This is 
because Br(r —y yrf) is proportional to tan 2 /3 while the box contribution does not depend on 
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tan/3. The double penguin contribution is proportional to tan 4 /3, and thus the amplitude can 
become larger for tan/3 > 30. For large tan/3 > 40, the constraint from Bt(B s — y /i/x) [38] 

Bi(B s ->■ fifi) < 4.3 x 10~ 8 , (19) 

becomes more important, since it is proportional to tan 6 /3. In the plots, the lines are terminated 
when the B s — y \i\i bound is saturated. 

The left-handed flavor changing Higgs coupling, X"[ 3 R can be generated by chargino diagram 
even if the flavor violating source is only the CKM mixing. Therefore, the mixing amplitude can 
be enhanced when only the right-handed down-type squark mass matrix has the off-diagonal 
element such as in the SU(5) case. The left-handed source of FCNC is also helpful to enhance 
the mixing amplitude since it can give a constructive contribution to the left-handed penguin. 
In fact, in the SO(10) boundary condition, there is additional phase freedom from the Yukawa 
matrix, and the phases of the off-diagonal elements for left- and right-handed squark mass 
matrices are independent in the basis where the down-type quark mass matrix is real and 
positive diagonal matrix. As a consequence, even if we fix the phase of M( 2 susy^ there remains 
one more phase freedom. In figure 2, actually, we choose the additional phase in the left- 
handed off-diagonal element to make the constructive contribution to the mixing amplitude. 
Therefore, the mixing amplitude under the SO(10) boundary condition can be larger than the 
case in SU(5). 

It is interesting to note that the chargino contribution to b — y s'j is destructive when the 
Higgs penguin contribution is constructive [151 [12]. This is roughly because the electric charge 
of down quark is negative, and the signatures of amplitudes for b — > S7 and the finite correction 
of the down-type quark mass matrix are opposite. As a result, the SO(10) boundary condition 
can be also preferable from the b — y 57 constraint. 

We comment on the GUT threshold effects for the boundary conditions, Eqs. (TT0]) . (TTT|) . In 
the SO(10) case, the flavor violation pattern in the lepton sector and the quark sector can 
depend on the SO(10) symmetry breaking vacua. Actually, in order to forbid a rapid proton 
decay, the quark flavor violation should be larger than the lepton flavor violation among the 
symmetry breaking vacua [39|. Namely, it is expected that Kq, ku c , an d Kd c are much larger 
than k l and %. For example, if only the Higgs fields (8, 2, ±1/2) are light compared to the 
breaking scale (which is the most suitable case), one obtains k,q = hjj c — k d c , and only quark 
flavor violation is generated, while the lepton flavor violation is not generated. On the other 
hand, when the flavor violation is generated from the minimal type of SU(5) vacua with the 
type I seesaw, it is expected that kl is always larger than k^c since the right-handed Majorana 
mass scale is less than the scale of colored Higgs mass. Therefore, the existence of b-s transition 
indicated by the experimental results in Fermilab predicts the sizable lepton flavor violation in 
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the minimal type of SU(5) model. In other words, if the results of a large B s -B s phase is really 
an evidence of NP, the minimal type of SU(5) GUT models are restricted severely [91 [TTJJ, ITT] . 

5 Constraints from the neutralino dark matter 

The cosmic microwave background anisotropy measurement by WMAP [TO] put a stringent con- 
straint on the SUSY parameter space through the dark matter requirement. Within 2 a, the 
neutralino relic density should be 0.106 < Qh 2 < 0.121. This assuming that dark matter con- 
sists solely of neutralino, i.e. smaller relic density cannot be a priori excluded. In SUSY models 
with universal gaugino and sfermion masses, m^j and m respectively, it is well known that 
there are five distinct regions that satisfy the relic density constraint: (a) the bulk region where 
both mi/2 and mo are small, (b) the neutralino-stau coannihilation region where the lightest 
stau mass is almost degenerate with the neutralino mass, (c) the focus point (FP) /hyperbolic 
region at large mo where the /i parameter becomes small and the lightest neutralino gets more 
Higgsino content, (d) the funnel region where the heavy Higgs masses (m^ and m#) is about 
twice the neutralino mass, and (e) the neutralino-stop coannihilation region where the lightest 
stop mass is suppressed by large off-diagonal terms when the trilinear coupling parameter Aq 
is large. In our analysis we assume that the Higgs soft masses are not tied to m , and from 
this assumption we have two more free parameters /i and m^. In such models there could be 
another dark matter region, in addition to the five regions above, i.e. the neutralino-sneutrino 
coannihilation region at large itla and/or /i [41] . 

The neutralino dark matter hypothesis is very attractive, and there are lots of activities, 
experimentally and theoretically, to discover the dark matter candidate. As a weakly interacting 
massive particle (WIMP), the lightest neutralino can in principle be detected directly by ultra- 
sensitive detectors. Such experiments are now reaching O(10 -8 pb) sensitivity level for certain 
values of neutralino mass JT2] IT3] . Neutralino particles in the galaxy can also give out some 
indirect signals through their annihilation, in particular from some regions where neutralinos 
can accumulate due to some gravitational potential attractors. It was pointed out that high 
energy neutrino flux from the sun can potentially be a clear signal of dark matter annihilation 
in the sun jTT], and this is currently being searched by the IceCube experiment and its upgrade 
with DeepCore which can lower the neutrino energy threshold for the detector [45J. 

To study the dark matter aspect of our model, we calculate the neutralino relic density, the 
neutralino-proton elastic scattering cross section, and the muon flux induced by solar neutrinos 
from neutralino annihilation. For the muon flux calculation we use the DarkSUSY program 
version 5.0.5 [46] which utilizes the results of WimpSim [47], interfaced with our own spectrum 
program. Solar neutrino from WIMP dark matter in various models has recently also been 
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Figure 4: (Left) The m^ — /i plane in SU(5) model with tan/3 = 40, m^j = 500 GeV, m = 
1 TeV, Aq = 0, and A^ p /A^ M = 1, showing various constraints from flavor and dark matter 
sectors as discussed in the text. (Right) Same plot for mi/ 2 = 800 GeV. 



analyzed in |IEJ SHI EQ]. We assume the NFW profile [51] for our analysis. As mentioned in 
[48J, there are uncertainties in the neutralino-nucleon cross section due to the strange quark 
role in the interaction. We use their default values for E^ an d A s , i.e. E^-jy = 64 MeV, and 
A® = -0.09. 

Since the parameter space in the minimal type of SU(5) is restrictive, (in other words, 
the mass spectrum is constrained), the discovery potential of this region appears to be very 
promising at the LHC, and at the direct and indirect dark matter search experiments |18| . Since 
small values of \x are not preferable due to the r — y fi'y constraint, the WMAP relic density 
prefers the funnel solution (i.e. the neutralinos annihilate through the heavy Higgs bosons pole). 
It is also true that a small value of m^ is preferred to enhance the double penguin contribution. 
In that case, the spin-independent neutralino-nucleon scattering cross section can be enhanced. 

In figure HJ we show the allowed region in the tua-^ plane when A^ F /Af^ = 1 in the SU(5) 
case, Ki = K£>c For the left panel we choose as SUSY parameters tan/3 = 40, m\/2 = 500 GeV, 
mo = 1 TeV and Aq = 0. The same parameters are used for the right panel except for 
rrii/2 = 800 GeV. The yellow and gray regions are excluded by the r — > /ry and B s —y jijj, 
constraints, respectively. The red-brown region is excluded because the lightest stau is the LSP 
there, and therefore the neutralino cannot be the dark matter candidate. The WMAP relic 
density range is obeyed for the neutralino in the narrow blue bands. The solid green lines are 
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contours for the scalar neutralino-proton elastic scattering cross section of 5, 1,0.1 x 10~ 8 pb 
respectively from left to right. The almost horizontal dashed green lines are for the spin- 
dependent cross section, 5, 10, 50 x 10~ 8 pb respectively from top to bottom. We also show 
contours of muon flux, labeled as Ex-y, where x is the assumed detector energy threshold in 
GeV and y is the flux in km~ 2 yr -1 . We show two cases for E(threshold): 100 and 10 GeV, 
dotted and solid respectively. As we can see, the Br(r — > /ry) constraint is important. In the 
left plot, the funnel happens at relatively small m^, and a large part is allowed. In the right 
plot, however, the funnel is shifted to the right due to the larger neutralino mass. In the later 
case, we need much larger // to satisfy the Br(r — > /17) constraint, although there would be an 
upper bound on fi due to the sneutrino LSP region jUJ. 

The muon flux rate is correlated to the neutralino-proton scattering cross section since this 
cross section determines the number of neutralinos accumulated in the core of the sun, hence the 
neutralino annihilation rate there. Since protons constitute a large portion of the sun, both the 
spin dependent and independent neutralino-proton cross sections are comparably important. 
For most of the MSSM parameter space the spin dependent part is larger than the scalar part, 
and this leads to a widespread misconception that for solar neutrino flux calculation only the 
spin dependent cross section is important. However, there are some regions of the parameter 
space where the scalar and the spin dependent cross section are about the same order of 
magnitude, and in this case the scalar contribution is also significant in determining the flux. 
This was also pointed out by [48J. We can see in Fig. H] that for small vha the scalar cross 
section is quite large, and the muon flux is following the scalar contour. For large tua, however, 
the spin dependent part becomes larger than the scalar part, and the muon flux contour is 
flatten out. Furthermore, the solar neutrino flux also depends on the branching fractions of 
the neutralino annihilation. Near the stau LSP region, the lighter stau mass is relatively small, 
enhancing the t + t~ channel, which in turn increases the flux (visible on the left panel). In the 
middle of the funnel region, the neutralino relic density is much below the lower bound of the 
WMAP range. We rescale down the muon flux due to this fact, and this is seen as a drops on 
the muon flux rate contour. 

The IceCube with DeepCore can potentially detect neutrinos with energy down to 10 GeV [52] 
It appears that for the left plot large allowed region, i.e., the entire left band of the funnel region 
and part of the right band, can be detected. Note, however, that we should also consider the 
backgrounds to have a clear discovery [49J , and therefore sufficiently large amount of data would 
need to be collected. For the plot on the right side, however, it would be very difficult for the 
IceCube with DeepCore to probe the WMAP region not already excluded by the Br(r — )■ /ry) 
constraint. 
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In the SO(10) case, the figures for the same parameter space remain unchanged qualitatively, 
except for a lower \x region where the chargino contribution of the Higgs penguin from the left- 
handed squark FCNC is important. As mentioned in the previous section, when the SO(10) 
breaking vacua is chosen to satisfy the proton decay constraint while gauge unifications are 
maintained, there is no stringent constraint from r — \ /17, and a larger region of parameter 
space is allowed. Consequently, the dark matter direct and indirect detection will play more 
significant roles in excluding the parameter space. 

6 Discussions 

In this paper, we investigated the effect of the recent dimuon CP asymmetry from B decay 
modes observed at 3.2 a deviation from the Standard Model by the DO collaboration in the 
context of -R-parity conserving SUSY GUT models and show that a large amount of flavor 
violation between the second and the third generation can be generated. Not only the large 
flavor violation arises due to large atmospheric mixing angle, but also new CP phases are 
obtained from the Yukawa interactions in grand unification, and they are responsible for the 
observed large CP asymmetry. 

Because of the quark-lepton unification, the CP asymmetry is restricted due to the r — > fi'-f 
bound. This restriction depends on the source of flavor violation (Dirac neutrino Yukawa cou- 
pling in the minimal type of SU(5) or Majorana neutrino Yukawa coupling in the minimal type 
of SO(10)), SUSY mass spectrum (larger SUSY breaking masses are preferred), and dominant 
diagram (box diagram for lower tan/3 or double Higgs penguin diagram for larger tan/3). We 
find that large values of tan /3 (= 30 — 50) are preferred because the CP asymmetry is enhanced 
via the double Higgs penguin diagram, whose contribution is proportional to tan 4 /3 and a siz- 
able contribution to the flavor violating b-s transition can be easily realized. We found that 
the minimal type of SO (10) is preferred due to the fact that both left- and right-handed squark 
mass matrices can have FCNC sources. The intermediate values of tan /3 (tan (3 = 20 — 30) are 
not very preferable. In the case of the minimal type of SU(5), tan (3 should be large to make the 
double penguin diagram dominant. In that case, B s — > \x\x is enhanced and it provides a lower 
bound of Br(B s —¥ /i/i), which is about 1 x 10~ 8 , in order to obtain a large CP asymmetry. 

The symmetry breaking vacua can be chosen to make the quark-lepton unification relaxed in 
the SO(10) case, while in the SU(5) case where the neutrino Dirac Yukawa coupling is the source 
of FCNCs, the leptonic FCNC is always larger than the quark FCNC. Therefore, the bound from 
the LFV in the SU(5) is more stringent, and in other words, the spectrum is more predictive to 
obtain the large CP asymmetry indicated by the like-sign dimuon charge asymmetry. In fact, in 
order to satisfy the dark matter content of the universe, the CP asymmetry prefers the funnel 
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solution where the lightest neutralinos annihilate through the heavier Higgs bosons pole. Such 
a restriction on the spectrum from the flavor physics would allow us to observe this parameter 
space at the LHC and at the direct and indirect dark matter detection facilities more easily. 
We showed that the high energy neutrino flux from the sun from the neutralino annihilation 
can be detectable at IceCube with DeepCore for some regions of parameter space. 
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